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ESTIMATING  MISSILE  RELIABILITY 
by 

S.  Blumenthal  and  J.  Denton 

I.  Introduction. 

The  problem  to  be  discussed  is  a  specialization  of  a  problem  men¬ 
tioned  by  Deemer  and  Mayberry  [l96l] .  Their  problem  concerned  the  allo¬ 
cation  to  targets  of  a  stockpile  of  missiles,  the  assignment  being  made 
on  the  basis  of  the  outcome  of  a  testing  program  designed  to  estimate 
missile  reliability.  Thus  the  question  of  hov  many  missiles  from  the 
stockpile  should  be  expended  in  operational  test  firings  must  be  studied. 
Somewhat  more  specifically  their  problem  may  be  put  into  the  following 
form.  We  assume  we  have  a  stockpile  of  missiles  and  a  set  of  targets 
^Tl’  T2’  **  *'  Tt^  over  which  we  wish  to  allocate  these  missiles  in  an 
optimal  way.  If  a  missile  is  completely  reliable  then  it  has  a  certain 

probability  P  of  destroying  a  target  T  at  which  it  is  aimed. 

J  J 

Since,  however,  missiles  are  not  completely  reliable,  we  assume  that  a 

missile  has  reliability  F  and  that  the  probability  that  the  target 

T.  survives  a  single  missile  is  (l-P.R),  while  the  probability  that 
J  J 

it  survives  n  missiles  is  (l-PjR)n.  However,  the  value  of  R  must 
be  estimated  frpm  the  operational  testing.  It  is  desirable  that  each 
target  upon  which  missiles  are  expended  should  have  enough  missiles 
allocated  to  it  so  that  its  survival  probability  is  very  small.  However, 
we  wish  to  avoid  assigning-more  missiles  than  necessary  to  a  given  target. 
Nevertheless,  a  given  amount  of  over-assignment  is  to  be  preferred  to 
the  same  amount  of  under-assignment;  i.e.,  in  terms  of  a  non-negative 


lose  function,  the  loss  is  zero  when  the  target  does  not  survive  and 
positive  otherwise. 

II.  The  Problem. 

The  specific  problem  to  be  discussed  below  may  be  formulated  as 
follows.  We  imagine  a  circle  of  fixed  radius  p  surrounding  each 
target  with  the  property  that  if  a  missile  detonates  within  this  circle 
the  target  will  be  destroyed,,  but  the  target  is  undamaged  by  a  missile 
detonation  outside  this  circle.  We  associate  with  each  missile  a  number 
R  ^  called  its  reliability  which  we  define  as  the  probability  that  the 
missile,  when  aimed  at  a  given  target,  will  detonate  within  the  circle 
of  radius  p  surrounding  that  target.  Thus  the  probability  that  a 
given  target  survives  one  missile  is 

1  -  R  , 

and  the  probability  that  it  survives  n  missiles  is 

(1  -  R)n  , 

so  the  probability  of  destroying  the  target  with  n  missiles  is 

1-  (1  -  R)n  . 

Let  us  confine  our  attention  to  the  target  T.  Given  a  number  Q, 
specified  in  advance,  with  0  <  Q  <  1,  the  number  of  missiles  n^,  to 
be  expended  on  the  target  T  is  then 


The  analysis  based  on  R  is  conceptually  the  same  as  if  R  is  replaced 
by  the  value  PR  developed  in  the  Introduction. 
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(1) 


r ! 


f' 


nT  «  mln(h:l^(l-R)h  >  Q)  . 


Since  R  is  unknown,  e  cannot  determine  nT.  We  therefore  wish 
A 

to  obtain  an  estimate  R  from  which  may  be  estimated  in  the  obvious 
way: 

(2)  n^  =  min{h:l*-(l-R)h  >  Q)  . 

The  loss  function  to  be  employed  in  this  situation  takes  the  value  zero 
if  T  is  destroyed  with  probability  greater  than  or  equal  to  Q  and 
one  if  this  probability  is  less  than  Q.  This  leads  to  the  search  for 

A 

an  estimator  R  which  underestimates  R  with  predetermined  probability. 
That  is,  given  a  ,  0  <  a  <  1,  we  seek  an  R  with  the  property  that 

p[r  <  R]  >  lwa. 

A 

In  classical  terminology  R  is  a  lower  confidence  bound  for  the  value 
of  R. 

We  also  assume  that  the  target  lies  at  the  origin  (0,0)  of  a 
Cartesian  coordinate  system.  If  . . .,  (xn,yn)  are  the  points 

of  impact  of  the  n  missiles,  then  each  is  thought  of  as  a  random  ob¬ 
servation  of  a  bivariate  normal  variable  with  mean  vector  (0,0)  (there 
is  no  aiming  bias)  and  covariance  matrix  £  given  by 


In  words,  X  and  Y  are  independently  normally  distributed  with  zero 

2  2 

means  and  unknown  variances  cr  and  cr^  ,  The  reliability  R  then 
has  the  form 
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(3) 


R  =  R(£)  =  PCX2  +  Y2  <  p2} 


where  p  is  a  given  positive  constant. 

We  now  consider  several  methods  for  obtaining  confidence  intervals 
for  R.  Once  these  are  obtained  we  employ  their  lower  bounds  as  values 

A  .  .  A 

R  to  be  substituted  into  equation  (2)  to  obtain  the  estimate  n^  . 

III.  Methods  of  Estimation. 

Method  1 

On  the  basis  of  the  sample  (x^,y^),...,  (xn>yn)  we  want  a  value 

A 

R  such  that  when  R  is  the  true  value  of  the  parameter,  the  probability 

a  , 

that  R  is  no  greater  than  R  should  be  at  least  (1-a): 

(4)  P{$  <  R|R)  >  1-a  . 

A 

The  experimenter  is  to  specify  a.  By  choosing  R  in  this  manner,  we 
are  essentially  "underestimating"  R  which  leads  to  a  conservative 
assignment  of  missiles. 

We  can  always  satisfy  (4)  by  setting  ^  =  0,  but  this  leads  to 
assigning  all  missiles  to  one  target.  Thus  we  should  like  to  employ  an 
estimator  ^  which  satisfies  (4)  while  maximizing  the  number  of  targets 
attacked.  In  classical  terminology,  this  means  that  we  are  looking  for 
a  shortest  upper  confidence  interval  for  R. 

Let  us  now  confine  our  attention  to  the  case  where  ai  =  a2  = 

In  this  case  we  have 
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(5) 


R  = 


2na 


a 


2  2  2 

-ifr  -£? 

e  dx  dy  =  1-e 


2/2.  2 
x  +y  <  p 


It  is  clear  that  R  is  a  monotone  decreasing  function  of  a,  so  that 

it  suffices  to  find  an  upper  confidence  bound  a  for  a  ,  i.e.,  the 

A 

desired  bound  R  is 


(6) 


A 

R  =  1 


2/2&£ 


Consider  the  following  hypothesis  testing  problem: 

Given  ^  =  \i2  =  0  ,  to  test 

„  2  2 
H0  :  °  -  °0 

2  2 

vs.  H.^  :  a  <  aQ  . 

It  is  well  known  that  (c.f.  Lehmann,  1959)  the  test  which  accepts 
n  (X^+Y  2) 

H.  when  2,  - 2 -  >  C]  is  a  uniformly  most  powerful  unbiased  test 

i=l  a 

of  Hq  vs.  with  size  a,  where  is  the  upper  lOO(l-a)^  point 

of  the  chi-squared  distribution  with  2n  degrees  of  freedom.  From  this, 
we  see  that  a  "uniformly  most  accurate"  upper  confidence  bound  'a2,  for 
a2  is  given  by 

(7) 


#  =  I  (X2  +  Y2)  . 
C1  i=l  1  1 


Note  that  the  bound  given  by  (7)  and  (6)  has  the  following  property. 
Among  all  bounds  $  satisfying  (4),  this  bound  minimizes  EpL(R,$) 
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where  L(R,6)  Is  any  measure  of  the  loss  resulting  from  overestimating 

R,  and  where  L(R,R)  is  non-negative  for  ^  >  R,  and  is  non-deereasing 

in  R  (see  Lehmann  [1959] >  p.  78).  is  then  obtained  from  (2). 

,  2  2 

Now  let  us  turn  our  attention  to  the  case  where  =  Ka^  ,  where 
K  is  a  known  positive  constant.  In  this  case, 


(8) 


'  2*0.^  /2  J2 


2  2 

+  -2_) 

'2  2 ' 

2  of  2 of 

e  dx  dy 


2^  2 
x  +y  <  p 


=  1 


"  "  f 

n  Jo 


n/2 


o*( l+(K-l)alne)de 


Again,  it  can  be  seen  that  R  is  a  monotone  decreasing  function  of 
and  again  an  upper  confidence  bound  ^  for  yields  the  lower 

confidence  bound  $,  where 


(?) 


2 

P _ 

_  p —  "  p" 

o*(l+(K-l)sin  e) 

1  d0  . 


This  last  integral  may  have  to  be  evaluated  by  numerical  methods. 

Proceeding  exactly  as  in  the  case  =•  a,  we  find  the 

a2  2 

following  most  accurate  upper  confidence  bound  for  : 


(10) 


^  o  o 

I  txf  +  <r  1  1 
A  1  K 


> 


where  is  the  upper  (l-a)$  point  of  the  chi-square  distribution 

with  2n  degrees  of  freedom. 
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A 

The  R  obtained  here  has  all  of  the  desirable  characteristics  which 
were  enjoyed  by  the  estimator  in  the  previous  case.  However,  it  remains 
to  assign  a  value  to  K.  This  might  be  done  on  the  basis  of  previous  ex¬ 
perience,  or  by  performing  preliminary  experiments  to  test  an  hypothesis 
of  the  form  H:  K  =  Kq  where  KQ  is  chosen  arbitrarily  (K^=  1,  say)  or 
on  the  basis  of  how  large  a  deviation  from  unity  is  significant. 

Extension  of  this  method  to  the  case  where  K  is-  unknown'  and  no  ■ 

assumptions  are  made  concerning  it  does  not  appear  to  be  feasible  since 

2  2 

we  may  then  formulate  no  hypothesis  corresponding  to  H:  a  >  for 
which  a  UMP  unbiased  test  exists.  A  similar  difficulty  intrudes  in  the 
case  of  non-zero  covariances. 


Method  2 


For  a  sequence  of  m  test  firings  let  the  random  variable  Xi  be 
one  or  zero  according  as  the  i  detonation  is  or  is  not  inside  the 
circle  of  radius  p  surrounding  the  target  (i=l,2, . . . ,m) .  Following 
Lehmann  (1959)  set 

m 

(11)  Y  =  l  X  +  U  , 

i=l 

where  U  is  independent  of  X^,  i=l,  ...,m,  and  has  a  uniform  distri¬ 
bution  on  (0,l).  (The  use  of  the  statistic  Y  is  equivalent  to  randomi¬ 
zation.)  Then  Y  has  probability  density 


(12) 


(i-fO-W 


0  <  y  <  m+1  , 


where  [y]  denotes  the  greatest  integer  less  than  or  equal  to  y.  The 

conditions  of  Lehmann's  Corollary  3  (Chapter  8  §5)  are  then  satisfied 
A 

and  R  is  then  the  solution  of 
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(13) 


PrR(Y  >  y]  =  a  , 


where  y  Is  the  observed  value  of  Y.  A  solution  exists  for  a  <  y  <  m+a. 

A  A 

For  m  +  a  <  y  we  take  R  =  1,  and  for  y  <  a  we  take  R  =  0.  This 
bound  is  then  uniformly  most  accurate  in  the  sense  of  Lehmann  and  has 
the  further  desirable  property  of  minimizing  E^lCr^R)  subject  to  the 
requirement 

(110  PrD [R  <  R]  >  1-a  for  all  R, 

For  large  samples  the  usual  normal  approximation  with  continuity  correc¬ 
tion  utilizing  the  statistic 

(15)  Y'  =  X  X 

i=l  1 

A 

may  be  employed.  Then  R  is  the  value  of  p  satisfying 

(16)  Pr[|  -  p)  >  a  /mp(l-p)  ]  =  1-a  , 

where  $  is  the  cdf  of  the  standard  normal  distribution  and  <t(a)  =  a. 
Method  3 

Since  it  is  not  possible  to  find  an  exact  test  of  the  hypothesis: 

Hq:  R  =  Rq  vs.  H^:  R  <  Rq  (whose  acceptance  region  provides  a  lower 

bound  on  R);  and  since  it  is  very  difficult  to  find  an  approximate 

test  of  this  hypothesis  (at  least  using  the  method  to  be  described)  we 

shall  describe  a  large  sample  test  of  H^:  R  =  R^  vs.  H^:  R  /  RQ.  From 

the  acceptance  region,  we  get  upper  and  lower  confidence  bounds  (^  j^u) 

A 

on  R  and  we  can  then  take  R  as  a  lower  confidence  bound.  If  a  is 

L 

the  significance  level  of  the  2-sided  test,  and  it  Is  symmetric  then 


! 

I 
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i 

! 
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1~2  should  be  the  confidence  level  of  the  interval  (fi^l).  Since  the 
test  we  use  is  not  UMP,  the  bound  obtained  will  not  be  sharp.  How  this 
bound  compares  with  those  obtained  by  the  other  methods  investigated  is 
not  known,  nor  is  the  actual  confidence  level  for  a  small  sample,  since 
asymptotic  distribution  results  are  used  in  deriving  the  test.  Now  if 
a  and  represent  upper  confidence  limits  of  coefficient 

(l-U-a)1/2),  then  R^,^)  is  an  a-level  bound  on  R.  This  is  the 

A 

crudest  and  simplest,  and  it  is  not  clear  that  the  bound  R^,  described 
below,  will  be  larger  than  R(  Cg) . 

2 

The  missile  landing  coordinates  X  and  Y  are  independent  N(0,  c^), 

p 

N(0, Og)  respectively.  We  want  a  bound  for 


R  = 


2*Va 


Li 


2  J2  2 

x+y  <  p 


dx  dy  . 


n  p  n  p  2  2 

I-et  u  =  £  \  ’  v  =  2  >  and  let  ai  *  ei  '  °2  =  V  ^ 

Max  Likelihood  Estimates  of  0.  and  are  -  and  ^  respectively, 

1  d  n  n 

where  u  and  v  are  sufficient  statistics  for  0^  ,  0g  and  have  the 

distribution  p(u,v,01,0g)  =  [Cn(uv)n'^2"1(01,0g)  n^2e  1^2(£-  +  ^-)]. 

Let  0°  and  0g  be  the  maximum  likelihood  estimates  of  under 

the  restriction  R  =  Rq.  Then  we  know  (see  Lehmann  (1959)/  PP«  510  -  311) 
that  under  Hq  (i.e.,  when  R  =  Rq)  that  (-2  log  A^)  (where  Ar  is  the 
"likelihood  ratio") 
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n  n 

.  /  N 2  -n  n/^O-Ox”  5  ,  1/U  .  v  \x 

An  =  (uv)  n  e  (S^)  exp{-  +  -^)) 

®1  02 

2  2 

has  for  large  n  approximately  a  distribution.  If  P{X^  <  )  =  1-T), 

then  we  would  accept  if  (-2  log  An)  <  C^,. 


Note  that  6^  and  0g  are  the  solutions  of 


(i) 

(ii) 


d 

d0. 


(p(u,v,e1,e2)  +  ^(©^©g))  =  0 


~  (p(u,v,e1,e2)  +  ur^^))  =  o 


and  R(0  ,0O)  =  R 


Looking  only  at  (i)  and  (ii),  we  see  that  0^,0g  can  be  found  as  func¬ 
tions  of  n,  so  (-2  log  An)  is  a  function  of  \x,  say  h(n).  We  can 

show  that  h(0)  =  0  and  h(-»)  =  h(+oo)  =■  ».  Thus,  if  h(n)  is  decreasing 

for  (i  <  0  and  increasing  for  (i  >  0,  there  will  be  exactly  two  values 
of  n  so  that  h(n)  =  CQ  ,  i.e.  h^)  =  h(n2)  =  Cq  (^  <  0  <  Hg)  , 
and  h(u)  <  C  for  m  <n  <  n_.  Therefore,  we  would  accept  H  for 

Cl  ±  c.  O 

any  ii  in  this  range.  Since  0^,0^,  the  solutions  of  (i)  and  (ii)  are 
functions  of  n,  R(©1,02)  is  a  function  H(n)  of  |i.  If  R  is  a 
monotone  (increasing)  function  of  n,  then  corresponds 

to  R(n^)  <  R( |i )  <  R(h2),  and  Hq  would  be  accepted  for  R(n  )  < 

R($  ,$2)  <  R(Hg)  where  0^0 2  are  the  solutions  of  (i)  and  (ii)  under 

the  restriction  R(01,0g)  =  Rq.  For  large  n,  this  test  has  size  a. 

It  follows  that  a  confidence  interval  of  coefficient  (l-,a)  is  given  by 
(R(*i. R(h2)),  and  a  lower  confidence  bound  by  R^). 
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The  above  reasoning  is  based  on  the  assumption  that  R(n)  is  mono¬ 
tone  increasing  and  that  h(ji)  has  the  properties  described.  Since 
^  h(n)  =  c  (i  R(n),  it  follows  that  h(n)  behaves  properly  if  and 
only  if  R(n)  is  monotone  increasing.  Thus  it  only  remains  to  show 
this  last  monotonicity.  It  appears  to  be  but  this  has  not  yet  been 
demonstrated . 

One  obvious  difficulty  about  comparing  this  bound  R(n1)  with 
other  possible  bounds  is  that  it  is  not  possible  to  solve  explicitly 
for  it,  and  a  numerical  solution  of  the  equations  (i),  (ii)  and 
h(n)  =  Cq  is  necessary. 

Method  4 

For  another  large  sample  procedure  for  obtaining  a  lower  confidence 
bound  on  R  when  the  X  and  Y  aiming  errors  are  uncorrelated  and  both 
means  are  zero,  we  may  adopt  the  following  approach.  Let 


.  1  £  +.  £\ 
2  2  2 

R  ..--A-  f  f  .  W 
0  a”V2  hk  2 

x  +y  <  p 


2  i  2  2  2  1  2  2 

51  =  ~  ^  Xi  '  s2  =  h 


S  -  1 


2  2  ’ 
x  y 

“2  +  2 

\81  S2 I 


2ltSlS2  J'2J2  2 

x  +y  <  p 


dx  dy 
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R  -4* 


1 

r  r  rVi 

_  8«a^a 

sl=al  12 

LL  I1  1 

x  +y  <  p 

2  2\ 

h  +  h) 

ai  °2 1 


dx  dy 


R2  “  >.  2 
os,. 


s2=a2 


8jta,  a: 


12  J2  J2  2 
x  +y  <  p 


2  . 
x  a  -21 

- e 


2  2 

2L.  +  JL 

2  2 

\Q1  a2 


dx  dy, 


Employing  the  constant  function  1  as  the  dominating  function,  a  straight¬ 


forward  application  of  Lebesgue's  dominated  convergence  theorem  shows  that 

2  _  ,  2  „  ,  ,  ..  „  dR  S& 

is^  ds|  Ss^s| 


is  continuous  in  and  b^.  The  continuity  of 


Os2)2 


and 


•s?A 

d  R 

2\2 

(os2; 


follows  similarly.  (The  justification  for  differ¬ 


entiating  under  the  integral  sign  is  found  in  a  theorem  on  page  67  of 
Math.  Meth.  of  Stat.  by  H.  Cramer).  Thus  $  is  asymptotically  normal 
since  the  conditions  of  the  following  theorem  of  Cramer  are  satisfied: 


Theorem  (Cramer  p.  366,  with  suitable  notation  changes).  If  in 

2  2  2  2  A 

some  neighborhood  of  the  point  s^  =  0^  ,  the  function  R  is 

continuous  and  has  continuous  derivatives  of  the  first  and  second  order 

2  2,  a 

with  respect  to  the  arguments  s^  and  s2  ,  the  random  variable  R  is 

asymptotically  normal,  the  mean  and  variance  of  the  limiting  distribu- 

2  2 
2  2  2  2 

tion  being  given  by  R  and  2(o^)  R^  +  2(a2)  Rg  =  v  respectively. 

Thus  /n  [~=-]  has  a  limiting  N(0,l)  distribution.  Now  set 


A 

V  = 


2(e2) 


dR  2,‘ 

-5  +  2(.a) 

5sl 
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1! 

]l 


VA  .A 

Sr  Sr  22 

We  know  that  and  — 3  are  continuous  in  s.  and  b_  and  that 

ds2  Ss|  1  2 


and 


a $. 


s rr>  a n  in  probability.  Therefore  — *  -»  R,  and  — »  -»  R0 

2  2  Ssf  1  *-2  2 


Ss! 


in  probability.  So  by  a  theorem  of  Slutsky  (p.  255,  Cramer)  v  -»  v  in 
probability.  We  now  have  recourse  to  one  more  pertinent  theorem: 


T1 

t 


V 

i  !  < 
i 


Theorem  (Cramer  §20.6,  p.  2^k).  Let  |  ...,  be  a  sequence 

of  random  variables,  with  the  d.f.s.  F^,  Fg, ...  .  Suppose  that  (F^(x) } 
tends  to  a  d.f .  F(x)  as  n  ->  »  . 

Let  q  ,qg, ...,  be  another  sequence  of  random  variables,  and 
suppose  that  {r^j  converges  in  probability  to  a  constant  c.  Put 


X  =  i  +  tj  ,  Y 
n  n  n  n 


^n  ’ 


Z  =  —  . 

n  ^n 


Then  the  d.f.  of  X  tends  to  F(x-c).  Further,  if  c  >  0,  the  d.f.  of 


n 


Y  tends  to  F(-),  while  the  d.f.  of  Z  tends  to  F(cx). 
n  c  *  n 


Now 


/s 

R-R 

v 


A 

R.  _  R 

A  A 

V  V 


R  R  / 

By  Slutsky's  theorem - >-  in  probability  as  n  -»».  By  Cramer's 

£  v  *  /V 

v  E  R 

§20.6  theorem  the  d.f.  of  -  tends  to  the  d.f.  of  —  .  So,  by  a 

A  .  V  7 

y  A 

R-R 

second  application  of  this  theorem  the  d.f.  of  — p-  tends  to  the 
a  a  v 

distribution  of  ,  i.e.  /n  has  a  limiting  N(0, l)  dis- 

v  v 

tribution.  Therefore  we  may  proceed  to  construct  a  large- sample  lower 
a- confidence  bound  as  follows 

rR-R 


Pr[  /n  (—7-)  <  b ]  =  $(b) 


A  $ 

PrfR  -  p  b  <  r],  =  *(b) 

Vn 
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where  *  is  the  standard,  normal  c.d.f.  and  b  Is  chosen  so  that 
♦(b)  =  1-a. 

We  remark  that  though  the  derivation  might  be  long  and  tedious, 
the  extension  of  this  method  to  the  case  of  non-zero  covariance  appears 
to  require  nothing  new  in  principle. 

IV.  Discussion. 

2  2  2  A 

When  =  Cg  =  cr  and  the  covariances  are  equal  then  if  R^  is 

the  estimator  of  R  given  by  Ml  (method  1}  and  Rg  the  estimator 

A 

of  R  given  by  M2,  we  have,  since  Rg  is  uniformly  most  accurate 
R'  <  R  =>  Pr(Rg  <  R'  }  <  Prjl^  <  R' )  . 

A 

We  note  that  Rg  enjoys  its  pleasant  properties  independent  of 
the  form  of  the  distribution  function  for  the  points  of  missile  impact 
(i.e.  whether  or  not  we  assume  the  coordinates  of  the  points  of  impact 
follow  a  bivariate  normal  distribution).  In  particular  this  is  true  in 
the  case  cr^  0g  when  a  bivariate  normal  distribution  is  assumed. 

Also,  for  large  samples  the  estimator  based  on  Y'  appears  to  involve 
less  computational  labor  than  those  introduced  in  M3  and  Mt. 

At  present,  because  of  the  approximations  involved  in  the  estima¬ 
tors,  it  seems  that  the  only  direct  way  of  comparing  M3  and  Mt  with 
each  other  as  well  as  with  Ml  and  M2  appears  to  be  by  use  of  a  Monte 
Carlo  method. 

With  respect  to  the  formulation  of  the  problem  employed  it  may  be 
remarked  that  the  form  of  the  loss  function  requires  assigning  enough 
missiles  on  each  target  attacked  to  insure  Its  destruction  before 


It 


including  additional  targets  among  those  to  be  attacked.  If  some  targets 
require  more  missiles  than  others,  then  with  a  limited  missile  stockpile 
those  requiring  fewest,  missiles  should  receive  first  assignment  and  so 
on  until  all  of  the  missiles  in  the  stockpile  are  assigned.  If  all  tar¬ 
gets  require  the  same  number  of  missiles  but  are  not  equally  important, 
this  may  be  reflected  by  assigning  higher  priority  targets  a  loss  greater 
than  1,  if  sufficient  missiles  are  not  assigned  to  insure  destruction 
with  probability  greater  than  or  equal  to  Q. 

It  might  be  noted  that  the  model  set  fouth  in  Section  II  may  be 
modified  in  several  ways  without  materially  changing  the  estimation 
problem  involved.  As  an  example  we  suggest  a  cumulative  damage  model 
in  which  the  target  is  surrounded  by  a  circle  CT  of  fixed  radius  p.  . 

Each  missile  causes  total  destruction  inside  a  circle  C._  or  radius 

M 

r  about  its  point  of  impact.  The  target  is  considered  to  the  totally 
destroyed  when  all  of  the  area  inside  the  circle  CT  is  destroyed  by 
missile  blasts,  i.e.  when  all  of  the  area  of  can  be  covered  by 

the  area  of  circles  of  radius  r  drawn  about  the  points  of  impact  of 
missiles  which  have  been  fired  at  the  target.  If,  once  again,  we  assume 
that  the  coordinates  of  the  points  of  impact  are  independently  distri¬ 
buted  and  have  a  bivariate  normal  distribution  about  the  point  of  aim 
2  2  2 

with  =  a2  =  a  ,  then  questions  about  the  number  of  missiles  to  be 

expended  to  obtain  a  given  percent  damage  with  a  certain  level  of  confi- 

2 

dence  may  be  answered  by  obtaining  an  estimate  for  o  .  We  might  al¬ 
ternatively  wish  to  fire  enough  missiles  so  that  the  expected  coverage 
of  the  circle  C^,  was  greater  than  some  specified  amount  (see 
Mbrganthaler  [1961]).  In  this  situation  it  might  also  be  well  to 


consider  whether  the  point  of  aim  should  coincide  with  the  target  for 
every  missile  in  order  to  obtain  maximum  coverage  from  the  missiles 
expended. 
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